Theory of transient spectroscopy 
of multiple quantum well structures 
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A theory of the transient spectroscopy of quantum well 
(QW) structures under a large applied bias is presented. An 
analytical model of the initial part of the transient current is 
proposed. The time constant of the transient current depends 
not only on the emission rate from the QWs, as is usually 
assumed, but also on the subsequent carrier transport across 
QWs. Numerical simulation was used to confirm the validity 
of the proposed model, and to study the transient current on 
a larger time scale. It is shown that the transient current is 
influenced by the nonuniform distribution of the electric field 
and related effects, which results in a step-like behavior of the 
current. A procedure of extraction of the QW emission time 
from the transient spectroscopy experiments is suggested. 



I. INTRODUCTION 

Transient spectroscopy of quantum-well (QW) struc- 
tures allows to study the emission processes from the 
QWs and thus to obtain information on QW parame- 
ters, such as the energy spectrum, phataionization cross- 
section, tunneling escape time, etc.pffl This technique 
is based on an analysis of the transient current or ca- 
pacitance relaxation upon the application of a large- 
signal bias across the QW structure. It complements 
admittance spectroscopy, which studies the ac current 
in the QW structure upon application of a small-signal 
voltage.B The transient spectroscopy of QWs has many 
similarities with the deep level transient spectroscopy 
(DLTS)H and enables a simple theory to be [derived for 
■ use in the processing of experimental data.H However, 
' there is an important difference in the carrier capture by 
deep levels and that involving QWs. In the latter case, 
the presence of a continuous energy spectrum for the in- 
plane motion in the QW allows the capture by emission 
of a single optical phonon rather Jjhan by a multi-phonon 
processes typical for deep levels.B As a result, the corre- 
sponding capture times are several orders of magnitude 
less than fot,deep levels and often do not exceed a few 
picoseconds!] This quantitative difference results in se- 
rious qualitative consequences. The processes of carrier 



transport between neighboring QWs can no longer be 
considered as infinitely fast. These processes may play 
a decisive role in the relaxation kinetics changing notice- 
ably the formulae of a simple theory, similar to the case 
of structures with very high concentration of deep levels. u 
In this work we present a theoretical description of the 
transient spectroscopy of QWs and discuss its possible 
applications. We obtain more general analytical expres- 
sions for the parameters of the transient current than 
those of Ref. [3[ The analytical model is confirmed by 
numerical simulation, and the procedure of the extrac- 
tion of the QW parameters from experimental data is 
discussed. 



II. ANALYTICAL MODEL 

We consider a QW structure containing M QWs (n- 
doped at sheet density Nd) of width L w separated by 
undoped barriers of width L b large enough to prevent 
inter- well tunneling (see Fig. 1). This structure is typical 
for quantum well infrared photodetectors (QWIPs).EI The 
QW structure is provided with a heavily doped (Ohmic) 
collector and a blocking emitter contact (for example, 
containing a Schottky barrier or p-n junction) which is 
often used to avoid DC current ancLthus to simplify the 
interpretation of experimental datau 

During the first period of the transient spectroscopy 
experiment, til the forward bias is applied to the emitter, 
and all QWs are filled by electrons with the equilibrium 
sheet density Nq rs Nd ■ In the second period, a large re- 
verse bias V is applied to extract electrons from the QWs, 
and transient current is recorded. The problem has some 
similarities with the treatment presented for the. kinetics 
of electron packet in a system of undoped QWs.E9. Imme- 
diately after the application of V, the electric field in the 
QW structure is uniform and given by E = {V + Vu)/L, 
where V&, is the built-in voltage between the emitter and 
collector, and L = ML W + (M + \)L b sa (M + \)L b is the 
structure thickness. This field causes fast removal of de- 
localized electrons from the structure at almost fixed Nq. 
This is a very fast_(on a ps time scale) component of the 
transient currenttiJ limited by carrier capture and transit 



times, which is manifested as an instantaneous current 
step in the case of limited time resolution of a measure- 
ment setup. We shall be interested here in a subsequent 
slow current relaxation caused by the QW recharging. 
The initial part of this relaxation (after completion of 
the fast transient) can be easily calculated. 

In the presence of external illumination, the emis- 
sion rate from the i-th QW with electron density Ni is 
GNi = a$Ni + 7exp[(e/ — e^/kT], where a is the pho- 
toionization cross-section, $ is the incident photon flux, 
£i is the QW ionization energy, et {Ni) is the Fermi energy 
of the electrons in the QW, and 7 is the thermoioniza- 
tion coefficient. In our estimates and further numerical 
simulations, we shall restrict our analytical calculations 
to the case of relatively low temperatures or high light 
intensities giving G = (7$. This restriction is not com- 
pulsory since all analytical formulae are applicable for an 
arbitrary relation between the optical and thermal gen- 
eration. 

We assume that the carriers emitted from the QWs 
drift with a constant velocity vj towards the collector. 
While traversing a QW, carriers are captured by the QW 
with a probability p (0 < p < 1). Carriers emitted from 
the fc-th QW give the following contribution to the carrier 
concentration in the i-th barrier (between i-th and (i+1)- 
th QWs): 



n« = — (l -*»)*-*, (i>k). 

Vd 

The total concentration in the i-th barrier n 



G 



^ ' 71J -<- ' 



i — k 



fc=l 



l-'d 



(1) 



(2) 



fe=i 



Since the change of a QW charge is determined by the 
balance between carrier capture and emission, we can, 
with the help of Eq. (0) , obtain the system determining 
the kinetics of all N: 



dN 

—^ = -GN t +n l _ 1 VdP = G 
at 



means that Vd is no longer constant but changes from 
point to point in an unknown way so that the behaviour 
of I(t) remains unknown. That is why we restrict our- 
selves to the initial stage of the slow relaxation when we 
can still assume that in the right-hand side of Eq. (0) 
N = Nq and Vd — const. This gives 
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and, hence, of the current in the external circuit I{i) 
which could be expressed in terms of N (t) '■ 
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In principle, we can obtain analytical (though rather 
cumbersome) solution of the linear system of Eq. ([|) for 
an arbitrary t. However, it would not be correct. The 
change in N causes re-distribution of the electric field 
and, hence, of the drift velocity in the system. This 
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where I e — eGNoM is the total emission rate from all 
QWs. Eqs. (||),(0) give us the relaxation time constant 
(inverse normalized slope of the current): 



( dl/dt \ - 1 _ P M - (1 - p)[l - (1 - p) M ] 



G[l-(l+pM)(l-p) M ] 



(8) 



In the most interesting case, when p<l and M ^> 1 
(which corresponds to practical QWIPs), parameters Iq 
and t are expressed as: 
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where g = \/(pM) is a transport parameter. If we char- 
acterize QW capture processes by the capture time r c 
or the capture velocity v c — L p /t c , which is.p£elated 
to capture probability as p — 1/(1 + «(i/t; c ),E3'tj then 
g = T c JT tr + 1/M w T c /r tr , where r tr = L/vd is the tran- 
sit time. Therefore, the parameter g corresponds to the 
photocurrent gain of a QWIPE 

It should be noted that the amplitude of the transient 
current Io is equal to the amplitude of the fast trapa 
sient (primary photocurrent) in a photoexcited QWIP.O 
In general, the time constant for the transient current r 
is determined not only by the emission time 1/G, but 
also by a transport parameter g, similarly to the case 
of DLTS for a very high concentration of deep centers. a 
Particularly, in the case j<C 1 we have r rs l/(gG) 3> 
1/G. Hence, one cannot obtain the photoionization cross- 
section from the transient spectroscopy experiment ignor- 
ing the correction factor (see Eqs. (|q)-(|lfj|)) dependent 
on g. Only in the limiting case g> 1 (or pM <C 1), the 
relaxation time tends to 1/G which corresponds to the 
simple model of Ref. 0. To fulfill this condition, one has 
to use QW structure with small capture probability and 
small number of QWs. In this case Io = 0.5 x I e , which 
corresponds to a high-frequency gain value of 0.5 for ex- 
trinsic photoconductors apd-QWIPs with large value of 
the low-frequency gain g.lHlEJIla Note that the capture 
probability p is a function of the electric field, decreas- 
ing with field, so that the simplified approach predicting 



t = 1/G can be accurate at high fields, but inaccurate 
at low fields. 

We point out that the value of the photocurrent gain 
g can be determined from the transient photocurrent 
in QWIP illuminated by a step-like infrared radiation, 
where the ratio of the amplitude of the fast transient 
to the steady^state photocurrent is equal to {1 — g[\ — 
exp(-l/g)]}li 



III. NUMERICAL SIMULATION 

The model presented above is justified only for the ini- 
tial part of the transient, since we neglected the mod- 
ulation of the electric field due to the QW recharg- 
ing. To check this model and to obtain a description 
of the transient current for a wider time interval, we also 
studied the transient processes using numerical simula- 
tion. A time dependent QWIP simulatoiOilj was used 
with a zero-current boundary condition for the reverse- 
biased emitter contact. We simulated the transient spec- 
troscopy experiment of Ref. || on GaAs/Al .25Ga .75As 
QW structure with the area S = 2 x 10~ 4 cm~ 2 con- 
taining 10 donor-doped QWs with L w =60 A QWs and 
Njj =5x 10 11 cm -2 , undoped barriers with Lt, =350 A, 
Schottky emitter contact (Vu = 0.75 V) and collector 
GaAs contact doped with donors at 10 18 cm~ 3 . The pho- 
toexcitatioruconditions were similar to those used in the 
experiment J3 

Figure 2(a) shows the transient current calculated for 
a reverse bias of 1 V, which for the given Vbt corresponds 
to the applied field Eq w40 kV/cm. The capture prob- 
ability was chosen to be p = 0.04 so that g = 2.5. The 
initial part of the transient (t < t' where r' ~ 43 /is is 
the position of the first step in Fig. 2(a)) is very well de- 
scribed by an exponential function with the amplitude 
and time constant calculated by our analytical model 
without any fitting parameters. The time constant ob- 
tained is r w 345 /is, while the emission time is somewhat 
smaller, 1/G = 300 /is. 

Starting from the time moment t = r', the transient 
current decays more rapidly, and displays a series of steps 
and shoulders. These features are due to the redistribu- 
tion of the electric field caused by the depleting of QWs 
(see Fig. 2(b,c)). The i-th step occurs when the electric 
field in the (M + 1 — i)-th barrier becomes zero. When 
this happens, the electron density in the (M + 1 — i)- 
th QW returns to its equilibrium value No, and this well 
does not contribute to the emission current. The electron 
transport in the region between this QW and collector is 
purely diffusive. Using Eq. rt|) and the condition of zero 
electric field in the Af-th barrier E = £^i ^f-jrp[ 
(seo is the dielectric constant), we obtain the following 
estimate for the time constant r'\ 
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For the case of Fig. 2 this estimate gives r' ~ 39 /is, 
which is in a good agreement with the results of numerical 
calculations (r' w 43 /is). 

Since t' ^ t (unless g is very small), only a small ini- 
tial part of the transient process is described by the expo- 
nential function with time constant r and amplitude Jo. 
Thus, the fitting of experimental data by an exponential 
function to extract the time constant r should be done 
over the interval < t < t'. The fitting over longer time 
intervals can result in a significant error in estimating r 
(the dashed line on Fig. 2(a) is an "intuitive" exponential 
fitting with the time constant r = 130 /is). It should be 
stressed that the measurement circuit should have RC- 
time constant (R is the load resistance and C is the QW 
structure high-frequency capacitance) much smaller than 
time constant r' for correct evaluation of r. 

To check the influence of the QW capture velocity v c on 
the time constant r, we simulated the transient response 
for different values of v c . The values of the amplitude and 
time constant of the initial part of the transient current 
extracted from numerical simulation are shown in Fig. 3. 
A good agreement between these results and formulas 
Eq. (||),(||) proves the validity of the analytical model. 

So far we have assumed that the photoionization cross- 
section (or emission rate) is independent of the local elec- 
tric field. While this is a good approximation for bound- 
to-continuum transitions, the photoionization cross sec- 
tion for bound-to-bound and bound-to-quasi-bound tran- 
sitions depends strongly on electric fieldB To investigate 
this effect, we compared the results of simulation for 
the cases of field-dependent and field-independent cross- 
section (see Fig. 4). We used the model of the field- 
dependent cross-section proposed in Ref. ||: 

a{E) = cr x 0.5crfc[(e 2 - eEL w /2)/(V2Ae)], (12) 

where cto is the optical cross-section, erfc(x) is the com- 
plementary error function, £2 is the ionization energy of 
the second level, and Ae is the variance of £2 due to fluc- 
tuationsll The values £2=4 meV and Ae = 3.5 meV 
also taken from Ref. f| were used. It is seen from Fig. 4 
that the field dependence of cross-section washes out 
steps on the I{t) curve. Moreover, the initial part of the 
curve, < t < t' , deviates significantly from an expo- 
nential function of the analytical model. Physically, this 
is caused by the decrease of the photoemission current 
from near-collector QWs due to the electric field redis- 
tribution. This effect makes the extraction procedure of 
the time constant r from the slope of the transient cur- 
rent more complicated. However, the transient current 
amplitude Iq is not affected by the field redistribution. 
Since the amplitude is directly related to the photoemis- 
sion current (Iq = 0.5/ e ) in the case g ~^> 1, we propose 
to use the amplitude of the transient current rather than 
its slope to extract the photoionization cross-section from 
experimental data. 



IV. CONCLUSIONS 



A theory of the transient spectroscopy of QW struc- 
tures is presented. Analytical expressions for the initial 
stage of relaxation current are derived. It is shown that 
the time constant of the transient current is a function of 
both the photoionization cross-section and the transport 
parameter g becoming ^-independent at g 3> 1 . Numeri- 
cal simulation is used to check the validity of the analyti- 
cal model and study the transient current in more detail. 
The procedure of extraction the QW emission rate from 
the experimental data is discussed. 

This work is supported in part by the NSF under Grant 
No. ECS-9809746. H. R. and A. S. also gratefully ac- 
knowledge the support of NSERC. 




FIG. 1. Schematic illustration of the conduction band dia- 
gram and the transient current in the transient spectroscopy 
experiment. 
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FIG 2. (a) Transient current, (b) electron density in QWs, 
and (c) electric field in the barriers for structure with 10 QWs. 
In (a), solid line - numerical simulation, dashed line - an- 
alytical model, and dotted line - exponential fitting of the 
transient current on a large time scale. In (b) and (c), labels 
indicate the index of the QW and barrier, respectively. 




Photocurrent gain g=1/(Np) 

FIG. 3. Dependence of the inverse normalized time con- 
stant G/t and amplitude of the transient current Io/I e on 
the value of gain g. Solid lines - analytical model, markers - 
numerical simulation. 
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FIG. 4. Comparison of the transient currents for the 
cases of the field-dependent (solid line) and field-independent 
(dashed line) photoionization cross-sections. The inset shows 
the field dependence of the photoionization cross-section. 



